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ABSTRACT 

We present an analysis of the spatial distribution of various stellar populations within 
the Large Magellanic Cloud. We combine mid-infrared selected young stellar objects, 
optically selected samples with mean ages between ~ 9 and ~ 1000 Myr, and existing 
stellar cluster catalogues to investigate how stellar structures form and evolve within 
the LMC. For the analysis we use Fractured Minimum Spanning Trees, the statisti- 
cal Q parameter, and the two-point correlation function. Restricting our analysis to 
young massive (OB) stars we confirm our results obtained for M33, namely that the 
luminosity function of the groups is well described by a power-law with index —2, and 
that there is no characteristic length-scale of star-forming regions. We find that stars 
in the LMC are born with a large amount of substructure, consistent with a 2D frac- 
tal distribution with dimension ~ 1.8 and evolve towards a uniform distribution on a 
timescale of ~ 175 Myr. This is comparable to the crossing time of the galaxy and we 
suggest that stellar structure, regardless of spatial scale, will be eliminated in a cross- 
ing time. This may explain the smooth distribution of stars in massive/dense young 
clusters in the Galaxy, while other, less massive, clusters still display large amounts of 
structure at similar ages. By comparing the stellar and star cluster distributions and 
evolving timescales, we show that infant mortality of clusters (or 'popping clusters') 
have a negligible influence on galactic structure. Finally, we quantify the influence 
of the elongation, differential extinction, and contamination of a population on the 
measured Q value. 
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1 INTRODUCTION 

This study is part a series to quantify the degree and evolu- 
tion of substructure within galaxies. Here we introduce the 
methods and use the LMC as a test case. The distribution 
of young stars in the LMC was first studied with objective 
algorithms by Feitzinger & Braunsfurth (1984), who found 
that the distribution was not random but had distinct pat- 
terns on all of the spatial scales probed (60 pc to 8 kpc). 
Other studies have come to similar conclusions, that young 
stars are distributed hierarchically within galaxies (e.g. Efre- 
mov 1984, Battinelli, Efremov, & Magnier 1996, Elmegreen 
& Salzer 1999, Elmegreen & Elmegreen 2003, Elmegreen et 
al. 2006, Bastian et al. 2005, 2007, Sanchez et al. 2007). This 
pattern observed in the young stellar distribution is thought 
to originate from the distribution and fragmentation of gas 
within galaxies (e.g. Elmegreen & Efremov 1996; Elmegreen 
& Falgarone 1996). 

It has been noted that many nearby star forming clus- 



ters also have hierarchical structure seemingly dictated by 
the structure of the dense gas of natal molecular clouds (e.g. 
Lada & Lada 2003). Using statistical techniques, Gutermuth 
et al. (2005) studied three Galactic clusters of varying de- 
grees of embeddedness and demonstrated that the least em- 
bedded cluster was also the least dense and the least sub- 
structured of the three. That result hinted at the idea that 
the youngest clusters are substructured, but that dynamical 
interactions and ejection of the structured gas contributes 
to the evolution and eventual erasure of that substructure 
in approximately the cluster formation timescale of a few 
Myr (Palla & Stahler 2000). By having an accurate model 
of the spatial evolution of stellar structures we can approach 
a number of fundamental questions about star-formation, in- 
cluding what is the percentage of stars born in "clusters" and 
whether this depends on environmental conditions. Using 
automated algorithms on infrared Spitzer surveys of star- 
forming sites within the Galaxy (e.g. Allen et al. 2007), such 
constraints are now becoming possible. 
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In addition to constraining models of star-formation on 
galaxy and cluster scales, the spatial distribution of stars can 
have important implications on a variety of other astrophys- 
ical problems. For example, the mixture of gas and stars can 
affect the temperature structure within Hn regions which in 
turn can influence metallicity determinations based on emis- 
sion line spectra (Ercolano, Bastian, Stasinska 2007). Thus, 
a small dependence on clustering properties of young stars 
on galactocentric distance (i.e. if there is a slight preference 
to form more concentrated stellar groups in the central parts 
of galaxies compared to the outer regions as may be expected 
if gas pressure plays a dominant role in cluster formation) 
could artificially enhance the metallicity gradient observed 
in spiral galaxies. 

Previous works on stellar structures have mainly con- 
centrated on the youngest stellar objects (OB stars), while 
the evolution of these objects have remained largely unstud- 
ied. Efremov (1995) studied the distribution of Cepheids, 
whose ages are significantly higher than OB stars, and 
found evidence of huge stellar complexes which he suggested 
may be made up of dissolved aged associations. Zaritsky 
et al. (2004) showed the stark contrast between the spa- 
tial distribution of young and old stellar populations in the 
LMC using two colour-magnitude cuts through the Magel- 
lanic Clouds Photometric Survey (MCPS). 

The goal of this paper is to quantify the timescale over 
which stellar structures within the LMC evolve. To this end, 
we use two methods to estimate the evolution and anal- 
yse each method's strengths and weaknesses. Our primary 
dataset comes from the MCPS where we make cuts in colour 
and magnitude space in order to select differing stellar pop- 
ulations (quantifying the age through comparison with stel- 
lar isochrones with a given initial mass function) which we 
compliment with catalogues of stellar clusters (Hunter et 
al. 2003) and embedded young stellar objects (Whitney et 
al. 2008). 

The paper is organised in the following way. In § [5] we 
introduce the observations and the methods used, while in 
§[3] we compare the properties (size and luminosity distribu- 
tions) of the young star-forming groups between the LMC 
and M33. In §[4] we extract twelve samples of stars of varying 
mean ages in order to compare their clustering properties, 
which is carried out in §[S] In §|§]we discuss the implications 
of our results and present a series of simple models to aid in 
the interpretation, and in §[7] we present our conclusions. 



2 OBSERVATIONS AND METHOD 

Our primary set of data comes from the Magellanic Cloud 
Photometry Survey of Zaritsky et al. (2004). This dataset 
consists of U, B, V, and I band imaging of essentially the 
entire optical extent of Large Magellanic Cloud. The original 
list consists of over 24 million sources, far too many for our 
algorithm, so we restrict our study to various cuts in colour 
and magnitude space. 

The method that we employ to first study the distri- 
bution of star-forming sites is the same as was presented in 
Bastian et al. (2007) for the study of the spatial distribu- 
tion of OB stars in M33. We refer the reader to that work 
for details and only provide a summary of the method here. 
Using the photometry of individual stars we make colour 



and magnitude cuts in order to select young stars (see § [3] 
for details of the cuts). We then construct a minimum span- 
ning tree (MST) of the spatial positions of stars which pass 
our criterifl The tree is then fractured (FMST) by apply- 
ing a breaking distance, Dbroak, in which all edges of the 
tree, longer than Dbroak, are broken (i.e. removed), creating 
independent groups. In order to avoid a selection bias which 
has hampered previous studies, we apply multiple breaking 
distances (20 in all) from pc to kpc scales. Throughout this 
paper, unless otherwise noted, we have used a value of five 
for Nmin, which is the minimum number of sources used to 
define a group. 

In § [4] we extend our study by not only looking at the 
youngest stars, but also the distribution of stars of various 
ages. 

Throughout this study we have assumed a distance to 
the LMC of 50 kpc (e.g. Marconi & Clementini 2005). For 
comparison with models (both fractal and uniform distri- 
butions) we will assume that the LMC is largely 2D - i.e. 
that it is a disk-like face on galaxy. This assumption will be 
justified in § 15.11 



3 SIZE AND LUMINOSITY DISTRIBUTIONS 

As a first step towards understanding the spatial distribu- 
tion of star-forming regions in the LMC we perform a sim- 
ilar study as was carried out on M33 (Bastian et al. 2007). 
We have used cuts in colour and magnitude space (namely 
B — V < 0.5 and Mv < —4.5 - uncorrected for extinction) in 
order to select young (massive) stars, resulting in 5093 stars 
which pass these criteria. The sample is then analysed using 
the Fractured Minimum Spanning Tree method, selecting 
20 breaking distances between the effective resolution limit 
of 1.5 pc and 1 kpc (with equal steps in logarithmic space 
of 0.1485). We refer the reader to Bastian et al. (2007) for 
details on the method. This method produces a catalogue of 
groups and their properties, namely size and luminosity, for 
each breaking distance applied. We also combined these cat- 
alogues into a 'total sample' from which we removed groups 
found by multiple breaking radii (i.e. duplicate groups that 
had the exact same position and radius) only retaining the 
groups found by the smallest breaking radius. We note that 
in the individual catalogues we did not remove duplicates 
(i.e. groups that were found by other Dbroak), as each break- 
ing distance was treated seperately. 

3.1 Size distribution 

We first look at the size distribution of the groups found. 
In Bastian et al. (2007) we found that for M33 there was 
not a characteristic size for star-forming regions, contrary 
to what is often quoted in the literature. We showed that 
previous claims of a characteristic size of ~ 100 pc for OB- 
associations (e.g. Bresolin et al. 1998) was due to the selec- 
tion of a single breaking distance or scale in which to carry 
out the analysis. 

1 An MST is formed by connecting all points (spatial positions 
in this case) in order to form a unified network, such that the 
total length (i.e. sum) of all of the connections, known as 'edges' 
or 'branches', is minimized, and no closed loops are formed. 



© 2008 RAS, MNRAS OOO.rjJfTTJ 



The evolution of structure in the LMC 3 



1000 




10 100 1000 

R (pc) 

Figure 1. The size (radius) distribution of the groups for three 
different breaking radii as well as for the total sample. Each distri- 
bution is well fit by a log-normal distribution whose peak depends 
on the scale probed. The dashed (red) lines show the best fitting 
log-normal function to each distribution. 
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Figure 2. The cumulative luminosity distribution of the total 
sample of groups. The solid line shows the best fit to the data of 
the form N(dL) oc L~ a dL, where the best fitting a and error are 
given in the panel. 



We carried out a similar analysis for the LMC, and come 
to similar conclusions. Figure [T] shows the cumulative size 
distribution of groups found with various breaking distances 
and for the total sample. The dashed (red) lines show the 
best fitting log-normal function to each distribution, which 
provides a good fit for each breaking distance. The total 
sample displays an excess tail at large radii (> 400 pc) rela- 
tive to the best fitting analytical fit, however it is well fit up 
until this point. This break is similar, although at slightly 
smaller scales (400 pc vs. 1 kpc) to that seen in the autocor- 
relation function of young stars in the LMC (Odekon 2008, 
her Fig. 3). It is also similar to the scale where a break in 
the HI distribution is observed, 180 — 290 pc (Padoan et 
al. 2001; Elmegreen et al. 2001; Kim & Park 2007). This 
feature, seen in the HI and in the young stars, has been 
suggested to relate to the scale height of the galaxy, where 
structures begin to change from largely 3D to 2D objects. We 
note that Odekon (2008) also found a break in the autocor- 
relation function for young stars in M33, whereas no similar 
scale was identified in the study of Bastian et al. (2007). 

We confirm the results found for M33, that the turn- 
over in the size-distribution depends on the breaking dis- 
tance used in an almost one-to-one way, suggesting that 
previous reports of a characteristic size in star-forming re- 
gions were due to the selection of a single breaking dis- 
tance. Hence, we confirm that the LMC does not display evi- 
dence for a characteristic scale in the star-formation process. 
Similar results were found for the spiral galaxy NGC 628 
(Elmegreen et al. 2006). 

3.2 Luminosity distribution 

The luminosity of each group was found by summing the 
V-band luminosity of all the stars found to be part of the 
group. For simplicity, we just concentrate on the 'total sam- 
ple' (i.e. all of the derived groups with multiple detections 
removed), although we note that the results for each break- 
ing distance were similar. The cumulative luminosity distri- 
bution is shown in Fig. [2] along with a power-law fit to the 
data of the form NdL oc L~ a dL, shown as a solid (red) line. 



We find that the distribution is well fit by such an analytic 
form, with a — 1.94 ± 0.03. This is in excellent agreement 
with what was found for the group luminosity function in 
M33 (Bastian et al. 2007). It is also in good agreement with 
the study of groups of OB stars in the SMC based on the 
MST method by Oey, King & Parker (2004). The implica- 
tions of this are discussed in detail in Bastian et al. (2007) . 
We note that the luminosity function may not directly relate 
to the mass function of the groups if there is a substantial 
age spread in the groups or if the initial mass function (IMF) 
of stars is dependent on environment (i.e. since we only sam- 
ple the upper end of the stellar IMF we are insensitive to 
any fluctuations below our sample limits). 



4 SPATIAL DISTRIBUTIONS OF 
SUBSAMPLES 

4.1 Age groups and galactic structure 

The above sections dealt solely with young OB stars, how- 
ever, using different colour and magnitude cuts in the CMDs 
of the LMC stars from the catalogue of Zaritsky et al. (2004) , 
we can select populations of different mean ages. This is 
done by taking colour-magnitude regions (corrected for fore- 
ground extinction of 0.41, 0.32, 0.25, 0.15 mag for U, B, V, 
and I, respectively) along the main sequence and red-giant 
branches of the CMD. In order to minimize the effect of 
larger contamination in the outskirts of the galaxy, and to 
assure complete geometrical coverage, we limited our selec- 
tion to the inner three degrees (in radius) of the LMC. 

For selection of the age boxes we follow the same 
procedure as Gieles, Bastian, & Ercolano (2008, hereafter 
GBE08). The adopted regions are shown in Figure [3] where 
the sizes of the boxes were chosen such that they enclose 
2000 stars within their boundaries. In order to determine the 
average age of each of the boxes, we constructed a synthetic 
population of stars, assuming a constant star-formation rate 
and a random sampling from a Salpeter initial mass function 
(1955). We used the Padova isochrones (Girardi et al. 2002 
and references therein, Z=0.008) to assign a magnitude and 



© 2008 RAS, MNRAS OOO.fTHTTl 



4 Bastian et al. 



Table 1. Colour /magnitude cuts and corresponding ages for 
different 'boxes' in the LMC (shown graphically in Fig. [3). 
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colour for each model star, and assume a distance to the star 
of 50 kpc. We then find the number of stars at each age which 
fall into our colour-magnitude regions, and take the mean 
age of the stars which are contained within the designated 
region. Two examples (boxes 2 and 6) are shown in Fig. [4] 
The selection cuts and mean ages are shown in Table [1] We 
note that using the median age of each box instead of the 
mean does not significantly change our conclusions. Since we 
are using a CMD to select stars, the primary parameter that 
we are selecting is in fact stellar mass and not stellar age. 
We implicitly assume that high and low mass stars (above 
a few solar masses) have the same spatial distribution. 

Once the boxes have been selected we can then look at 
the spatial positions of the stars for each box. Six examples 
are shown in Fig. [5]where the box number and mean age are 
given in each panel. It is clear that the younger populations 
contain considerably more substructure than the older pop- 
ulations. The LMC bar is clearly visable in the older boxes, 
namely boxes 7, 9, and 12 in Fig. [5] and is clearly dominat- 
ing the spatial distribution in box 12. In the next section we 
will quantify these differences. 




SPATIAL ANALYSIS 
STRUCTURES 



QUANTIFYING 



Figure 3. Colour-magnitude diagram used to select the boxes. 
The size of the boxes was determined by requiring 2000 sources 
in each colour-magnitude box. Stellar isochrones from the Padova 
group (Girardi et al. 2002) are shown for logarithmic stellar ages 
(in years) of 7, 7.5, 8, 9, & 10 and Z=0.08 metallicity. 



5.1 Q parameter 

Cartwright & Whitworth (2004) have developed a formal- 
ism to quantify the structure of young star clusters. Their 
method uses the normalised mean MST branch length, m, 
and the normalised distance between all sources within a re- 
gion, s. The ratio between these two quantities, m/s= Q, is 
able to distinguish between a power-law (centrally concen- 
trated) profile and a profile with sub-structure (i.e. a fractal 
distribution). Additionally, this parameter can quantify the 
index of the power-law or the degree of sub-substructure, 
which, if one assumes is due to a fractal nature, its fractal 
dimension can also be estimated. Assuming a three dimen- 
sional structure, if Q is less than 0.79 then the region is 
fractal, larger than 0.79 refers to a power-law structure, and 
a value of 0.79 implies a random distribution (a 3D fractal 
of dimension 3 is a random distribution). However, if the 
distribution is two dimensional, as is approximately true if 
one is looking at a disk- like galaxy face on (like the LMC), 
then a random distribution has a Q value of 0.72. 

While this method was developed and has been applied 
to young embedded clusters (Cartwright & Whitworth 2004, 
Schmeja & Klessen 2006, Kumar & Schmeja 2007), it has 
also recently been applied to full stellar populations within 
the SMC (GBE08), finding that the Q parameter increases 
with the mean stellar population age. 

During our analysis, we found an unexpected relation 
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Figure 5. The spatial positions of the stars in boxes 1, 3, 5, 7, 9, & 12. Only the inner three degrees in radius are used in the analysis. 
The bar of the LMC is clearly visable in older populations, and dominates in age box 12. 



between the measured Q value and the elongation parameter 
introduced by Schmeja & Klessen (2006), in that more elon- 
gated clusters have a lower Q-value. We note here, that while 
this does not affect our results (each box has an elongation 
less than 1.2) it may be a serious problem for young embed- 
ded clusters which often have a filamentary structure (e.g. 
Gutermuth et al. 2008). This issue and a proposed solution 
are dealt with in more detail in Appendix A. Additionally, in 
Appendix B we discuss the effects of differential extinction 
on the Q parameter. 

In the bottom panel of Fig.[6]we show the measured Q- 
value for each of our 12 boxes (filled circles) . The Q-value is 
clearly increasing with age, beginning at a value of Q = 0.57. 
In order to measure the age at which the spatial distribution 
becomes uniform, we fit the Q-value as a function of the log- 
arithm of age (shown as a dashed line in the bottom panel 
of Fig. [6] and calculate when this function crosses that of a 
flat distribution (i.e. when the Q- value becomes constant). 
This happens at ~ 168 ± 30 Myr, where the error was cal- 
culated by the uncertainty in the value of the flat region of 
the diagram, with a value of 0.73. 

This is similar to what is expected for a 2D slightly cen- 
trally concentrated stellar distribution (a 2D power-law dis- 
tribution with index 0.3, decreasing outward; Cartwright & 
Whitworth 2004) , as is our reference distribution (see § 15. 2|) . 
This indicates that the spatial distribution of the stars be- 
comes more homogeneous as the populations age. However, 
box 12 has a much larger Q-value, which we attribute to the 
LMC bar, a strongly centrally concentrated 3D structure 
in the galaxy. Errors were estimated by generating artificial 
populations with the same number of sources, deriving their 
Q-value, and taking the standard deviation. 



5.2 Two point correlation function 

The second technique that we used to determine the 
timescale for structures to be erased in the LMC is the two- 
point correlation function (TPCF) . This method determines 
the distance between all possible pairs of stars (d s ), shown as 
a histogram, which is then normalized to that of a reference 
distribution, i.e. A r n n k s (d s ) / A rcfcronco (d s ). Here, A rcfcronco 
is taken as taken to be a slightly centrally concentrated pow- 
erlaw with index 0.3, which is a good representation of the 
inner three degrees of the LMC for stars with V < 20. This 
is similar to that done by Gomez et al. (1993) who used the 
TPCF to study the distribution of pre-main sequence stars 
in Taurus, although these authors used a random distribu- 
tion as their reference distribution. We are then left with 
a histogram for each age box of the number of connections 
with a certain distance, relative to the reference distribu- 
tion. For each of the age boxes we fit the distribution (in a 
logarithm of distance vs. logarithm of the number) with a 
linear relation to get the slope and zero-point. 

If the distributions are becoming more similar to the 
reference distribution, then the slope and zero-points should 
tend towards values of (see Appendix [C] for a demonstra- 
tion of this). The results for all age boxes are shown in the 
top and middle panels of Fig. [6] Similar to the Q-parameter 
method, we find that the distribution for the boxes is ap- 
proaching the reference distribution as the mean stellar age 
increases, for both the zero-point and slope of the distribu- 
tions. The fact that the values do not go precisely to zero 
is most plausibly due to a slight mis-characterization of the 
reference distribution. Errors on each of the points were cal- 
culated by randomly selecting 1100 stars in the age box and 
repeating the analysis. 

As can be seen by comparing all three panels in Fig. [6] 
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Figure 4. Examples of the age dating process. The histogram 
shows the fraction of stars which fall in our colour-magnitude 
selection box. The dashed line shows the mean age of all the stars 
which fall in the colour-magnitude box. See text for details. 

the three measurements of the disappearance of substructure 
agree very well. For the two measurements provided by the 
TPCF, we fit the zero-point and slope as a function of the 
logarithm of age (shown as dashed lines), and find that this 
crosses the flat portion of the diagrams at 175 ± 25 Myr, 
where the error indicates the uncertainty in defining the flat 
region of the diagram. 

5.3 Comparison of the methods 

The three measurements^ of the erasure of substructure 
within the galaxy (the Q-parameter, and the slope and zero- 
point of the TPCF) are all in close agreement. 

These results show that stars are formed with a large 
amount of substructure in the LMC and that, as this popula- 
tion evolves, it becomes more uniform, becoming statistically 
indistinguishable from the background reference distribution 
within ~ 175 Myr. Through the Q-parameter, and compar- 
ing with artificially generated distributions, as was done in 
Cartwright & Whitworth (2004), we estimate that stars are 
born with a fractal distribution, with a 2D fractal dimension 
of ~ 1.8. 

2 While we have three measurements, we caution that only two 
are independent, as the slope and zero-point of the TPCF are 
intimately related. 



Figure 6. The results of the three measurements (two techniques) 
of the erasure of substructure within the LMC. The top and mid- 
dle panels show the measured zero points and slopes of the two 
point correlation function for all age boxes. As the reference dis- 
tribution we use a centrally concentrated power-law distribution 
with index 0.3. Note the smooth transition from younger (more 
substructure) boxes to older ones. Both the zero point and slope 
distributions become flat at ~ 175 Myr, suggesting that most 
substructure has been erased by this age. The results of the Q- 
parameter are shown in the bottom panel, which largely agree 
with the results of the TPCF, i.e. that structure is erased in 
~ 170 Myr. Additionally, we show the evolution of Q for the 
cluster sample of Hunter et al. (2003). 



Such evolution of stellar structures has also been seen 
in NGC 1313 (Pellerin et al. 2007) and NGC 4449 (Anni- 
bali et al. 2007) using HST-ACS data. Both studies found 
that massive stars (OB stars) contain more structure than 
older stellar populations. Pellerin et al. (2007) explain this 
evolution of structure as due to 'infant mortality', i.e. that 
young stars are often formed in dense clusters which then 
'pop' (Kroupa 2002) when the residual gas left over from 
the star-formation process leaves the cluster (e.g. Bastian & 
Goodwin 2006). Due to our ability to quantify the timescale 
and degree of clustering we can rule out this scenario, which 
we will return to in § 16.11 

In such statistical analyses as presented here, contami- 
nating sources may potentially seriously affect the results. In 
order to estimate the magnitude of these effects we have run 
a series of simulations where we generate artificial distribu- 
tions with known parameters, and substitute random (or ref- 
erence) distributions with a given percentage. The results of 
these experiments are shown and discussed in Appendix [Cl 
In particular, we note that young stars will contribute to 
all age boxes (with the exception of box 12), and since they 
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display a larger degree of substructure, some degree of sub- 
structure is expected to be present in all boxes. However, 
the contamination is expected to be low (much less than 
30%), hence the statistical methods should not be affected 
by them, although they can be seen by eye in Fig. [5] 



6 COMPARISON OF OPTICAL GROUPS, 
YOUNG STELLAR OBJECTS AND 
STELLAR CLUSTERS 

In order to further understand the evolution of stellar struc- 
tures within the LMC we take advantage of two recent sur- 
veys. The first is that of Whitney et al. (2008) who used the 
SAGE (Meixner et al. 2006) multi-band Spitzer survey of the 
LMC to select high probability young stellar objects (YSOs) . 
These populations are expected to trace the evolution of 
forming stars/clusters from the parent clouds, through the 
embedded phase to the point where they can be detected at 
optical wavelengths. The second survey was that of Hunter 
et al. (2003) who catalogued the cluster population of the 
LMC, and measured spatial positions, magnitudes, and ages 
for over 900 clusters. 

Constraining the YSO sample to the inner three de- 
grees of the LMC (as was done for the optically selected 
samples) we find Qyso = 0.68. This is significantly higher 
than that found for the young optically selected (OB star) 
sample (Qob = 0.57). Since AGB stars and YSOs have sim- 
ilar colours even in the mid-IR (e.g. Whitney et al. 2008) it 
is possible that this difference is due to contamination of the 
YSO sample. If it is due to contamination (see Appendix [Cj) 
then this implies that the YSO sample has approximately 
40% contamination from stars which have a similar distri- 
bution to that of our reference distribution. This may be 
the case if AGB stars are the main contamination. The al- 
ternative interpretation, is that the YSO sample is made up 
largely of lower-mass stars than our OB star sample (Whit- 
ney et al. 2008), and that the distribution of star formation is 
mass dependent. Such extreme spatial dependence on mass 
appears to be at odds with observations of young Galactic 
star-forming regions (as indicated by the similarity of the 
stellar IMF in different environments) . Finally, we note that 
the lower spatial resolution of the YSO sample may act as 
to raise the measured Q-value. 

For the cluster distribution, we again limit the sample 
to the inner three degrees in order to draw useful compar- 
isons with our optically selected samples, and additionally 
provide a magnitude cut in the data of My < —4.0 in order 
to avoid incompleteness effects in the central region. We then 
break the cluster sample up into 8 age bins, with each bin 
containing an equal number of clusters, namely 69. The re- 
sulting distribution of Q-parameters vs. cluster age is shown 
in the bottom panel of Fig. HJ] As seen in the stellar distribu- 
tions, the Q-parameter of the cluster populations increases 
with mean age, and becomes flat between 100 — 200 Myr. 
The larger scatter in the cluster points relative to the stellar 
boxes is due primarily to the lower number of points used 
to calculate Q (69 vs. 2000). Overall, the cluster population 
loses structure in a similar way to the stellar population, and 
if anything, the structure removal is faster for the clusters. 
We will return to this point in § 16.11 



6.1 Interpretations: Infant mortality vs. general 
galactic dynamics 

6.1.1 Infant Mortality 

It has been shown that the majority of clusters that form be- 
come disrupted within the first few 10s of Myr of their lives 
(e.g. Lada & Lada 2003). The mechanism believed to be re- 
sponsible for this is the rapid removal of gas, due to stellar 
winds and/or supernovae, from the natal cluster. This gas 
removal leaves the cluster in an unbound state due to the 
removal of a large fraction of the gravitational potential of 
the cluster, that leaves the remaining stars in a super-virial 
state. If the star-formation efficiency is low (< 30%) the en- 
tire cluster may become unbound and expand rapidly. These 
types of clusters have been referred to as 'popping clusters' 
(Kroupa 2002) and evidence of their expanding nature has 
been seen in their low light-to-mass ratios for a given age 
(Goodwin & Bastian 2006). The clusters are predicted to 
expand with a velocity close to their velocity dispersion and 
may have a severe effect on the observed galactic structure. 

The standard scenario for the 'popping' clusters is that 
the stars and gas are in, or close to, virial equilibrium when 
the natal gas is removed. Thus, we can use the virial equation 
to find the relation between a cluster's mass (gas+stars) and 
velocity dispersion (e.g. Spitzer 1987). Due to the relatively 
low star-formation rate in the LMC we do not expect a large 
population of massive clusters. Assuming a stochastically 
sampled power-law mass function of the clusters (gas+stars) 
with index —2, lower mass limits of Mi ow = 1, 10, & 100 Mq, 
and cluster radii of 1 pc, we can set up a simple model of the 
evolution of the structure of a galaxy due to infant mortality. 

We start with a spatial distribution of clusters with the 
same Q-value as the youngest observed clusters. At each of 
these cluster centres we place a cluster of mass M c i and ra- 
dius 1 pc We then assign a number of stars to each cluster, 
directly proportional to its mass, with a gaussian spatial 
distribution (stochastically sampled). The velocity disper- 
sion of each cluster (gas+stars) is derived through the virial 
equation, assuming a given star-formation efficiency (SFE, 
c.f. Goodwin & Bastian 2006) . The cluster expansion is mim- 
icked by increasing R c i, with R c i= tdi s 0"ci, where fdi s is the 
time since the cluster began expanding and <r c i is the cluster 
velocity dispersion. At each timestep, we select 1000 stars at 
random and calculate their Q-value. Only stars within the 
inner 3 degrees (the simulations were set up relative to the 
observed LMC distribution) were selected. 

The results of these simulations are shown in Fig. [7] 
where the top, middle and bottom panels represent simu- 
lations where the mass function of the clusters is sampled 
down to 100, 10, & 1 Mq respectively. Since we have used 
the same number of clusters in each simulation, the horizon- 
tal shift in the simulations is related to the average velocity 
dispersion of the clusters, which decreases with decreasing 
Mi ow . If clusters start more compact (as suggested by obser- 
vations of young clusters, e.g. Mackey et al. 2007, Bastian 
et al. 2008), this would shift the curves to the left. 

The cluster mass function is thought to be continuous 
down to single stars (e.g. Oey et al. 2004), arguing for the 
applicability of the simulations with Mi ow = 1 Mq. This 
suggests that the erasure of structures is not being driven 
by 'popping clusters'. However, a caveat to these simula- 
tions is that we have assumed that all clusters pop in the 
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same manner, regardless of mass. This, however, may not 
be true, as the number of O-stars (which drive the gas re- 
moval, and hence the 'popping') in the cluster may influence 
the expansion properties (e.g. Baumgardt & Kroupa 2008). 
Additionally, we have not taken the motion of the clusters 
through the galaxy into account which also destroys sub- 
structure (see § 16. 1 . 2 f) . 

We can test the effect of infant mortality another way, 
namely by looking at the evolution of the Q-value for stars 
and that for catalogued star clusters. If infant mortality was 
a dominant cause of structural evolution within the galaxy, 
then we would expect the Q-values of the stars to increase 
more rapidly with age than that of clusters. This is because 
the star clusters which survive the 'infant mortality' phase 
will not have their positions changed, whereas the stars that 
formed in clusters that 'popped' will be expanding away 
from their natal positions. In other words, the stars will 
have two components to their velocity, namely expansion 
and motions through the galaxy, whereas clusters will only 
have one component, the velocity dispersion of the galaxy. 

For the cluster catalogue we use the Hunter et al. (2003) 
survey presented in §[6] In the bottom panel (open circles) of 
Fig. [6] we show the measured Q-value for each of the cluster 
sub-samples. In general we see the same behavior as found 
for the stellar boxes, namely that the Q value increases with 
mean age, beginning with a highly fractal, or substructured 
distribution and becoming flat near the same value as that 
found for the stellar boxes. The evolution between the stars 
and clusters is quite similar, in fact the clusters seem to 
evolve slightly faster (although this does not appear to be 
statistically significant given the small number of clusters 
per age bin) than the stars. Hence we conclude that the same 
mechanism is erasing the substructure seen in the stellar 
and cluster distributions, which eliminates the possibility 
that 'popping clusters' are the dominant cause, although 
this effect is likely taking place within the galaxy. 

6.1.2 Galactic Dynamics 

We conclude that the dominating effect of structure evo- 
lution in the LMC is general galactic dynamics. It is the 
velocity dispersion of the galaxy which eliminates substruc- 
tures and pushes the stellar distributions toward a random, 
or background, distribution. 

These results are consistent with the findings of the dis- 
tribution of Cepheids in the Galaxy (Efremov 1995) and the 
LMC (Marconi et al. 2006). For the Galaxy, Efremov (1995) 
found evidence for the presence of large 'Cepheid clouds' 
with sizes of a few hundred parsecs, which, considering the 
age of Cepheids (a few 10s to 100s of Myr) is what is ex- 
pected from our results. More directly, Marconi et al. (2004, 
their Figs. 4 & 5) found that young Cepheids (20 - 80 Myr) 
display much more substructure in their spatial distribution 
than older Cepheids (> 80 Myr) whose distribution is be- 
coming closer to random, in agreement to that found by 
Elmegreen & Efremov (1996). Finally, we note that this is a 
similar timescale (~ 150 Myr) to that found by Efremov & 
Elmegreen (1998) for the average separation between pairs 
of star clusters in the LMC to become constant, i.e. for the 
cluster population to become statistically uniform. 

The LMC has a velocity dispersion of ~ 22 km/s (Car- 
rera et al. 2008), if we adopt a radius of the LMC of 



3 — 4 kpc this translates to a crossing time (t C ross = R/&) of 
~ 135 — 180 Myr. This is remarkably similar to the timescale 
derived for the elimination of structure within the LMC. 
From this, we conjecture that structure will be erased on a 
crossing-time, independent of spatial scale. This may explain 
why young lower-mass clusters often contain large amounts 
of substructure (e.g. Gutermuth et al. 2005), while many 
of the massive young clusters (such as Orion) have much 
smoother central regions, where the crossing time is signifi- 
cantly lower than the outer regions (see Elmegreen 2008). 

In GBE08 we show the same quantitative relationship 
between the stars and clusters in the SMC and reach a simi- 
lar conclusion, namely that popping clusters are not driving 
the evolution of substructure within the galaxy. In the SMC, 
substructure is removed in ~ 80 Myr, which is again similar 
to the galactic crossing time. Note that the erasure of struc- 
ture in the 'popping cluster' scenario (e.g. Fig. [7} is mainly 
dependent on the average internal velocity dispersion a" me an 
of the clusters. Since <7 mca n is independent of the number 
of clusters in a population (it only depends on the mass 
function index and lower mass limit of clusters - which we 
assume to be invariant) then all galaxies should have struc- 
ture erases in approximately the same time. However, in the 
timescale of structure erasure is twice as fast as in the LMC. 

An interesting difference between the SMC and LMC 
stellar and cluster distributions is the Q-values which they 
attain when the become flat. For the SMC the Q-value of 
both clusters and stars becomes flat at Q = 0.81 whereas 
for the LMC, both distributions become flat at Q = 0.73. 
This difference shows the ability of the method to quan- 
tify the structure of distributions, as the SMC is largely 
three dimensional whereas the LMC is mostly disk-like, or 
two dimensional (both slightly centrally concentrated). The 
expected Q-values for homogeneous 2D and 3D distribu- 
tions are 0.71 and 0.79, respectively (Cartwright & Whit- 
worth 2004), therefore explaining the observed difference 
between the Q-values in the LMC and SMC. 



7 CONCLUSIONS 

We have presented an analysis of the spatial distribution 
of stellar structures within the LMC. By using the Zarit- 
sky et al. (2004) Magellanic Cloud Photometric Survey and 
sampling different parts of the stellar colour-magnitude dia- 
gram, we have selected stellar populations of varying mean 
age. This is supplemented with other datasets, namely the 
Whitney et al. (2008) Young Stellar Object (YSO) infrared 
sample and the star cluster catalogue of Hunter et al. (2003) . 

Based on the analysis of the distribution of OB stars, 
we confirm our earlier results on M33, namely that the star- 
forming groups are well approximated by a power-law lumi- 
nosity function with index of —2 and that there is no pre- 
ferred size in the distribution of star-forming sites (e.g. that 
of the OB association) . Instead, we find a continuous distri- 
bution of sizes, and find that previous reports of a character- 
istic size were due to resolution effects, along with choosing, 
apriori, a single size-scale to probe. 

We have used two methods to quantify the amount of 
substructure within a galaxy, the Q method (Carthwright 
& Whitworth 2004), and the two-point correlation function 
(TPCF). We found that stars are formed in a highly struc- 
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Figure 7. Simulated evolution of a galaxy due to infant mortality 
for various star formation efficiencies (SFE), see text for details of 
the simulations. Three values of the lower mass limit of the cluster 
mass function are given. All simulations begin with clusters with 
radii of 1 pc. The observed change in Q vs. time is shown as filled 
circles. 



tured way (approximating the distribution as a 2D fractal we 
find the fractal dimension to be ~ 1.8) and that this distri- 
bution evolves rapidly, becoming indistinguishable from the 
background or reference distribution after ~ 175 ± 30 Myr. 
We also found that the distribution of stellar clusters follows 
that of the stars of similar age, adding confidence in our age 
dating method as well as our structural analysis. 

Thus, our findings suggest that stars are formed in a 
highly substructured distribution and that this distribu- 
tion evolves rapidly, with small scale structures dispersing 
first, and that all of the original structure is removed in 
~ 175 Myr. This is approximately the crossing time of the 
galaxy and as such, we suggest that structures in galaxies, 
presumably regardless of spatial scale, are eliminated in a 
crossing time. Similar analysis done on the SMC (GBE08) 
finds that substructure is erased on ~ 75 Myr, which is com- 
parable to the crossing time of that galaxy. 

By comparing the observations to simple models, we 
conclude that infant mortality (or 'popping star clusters') 
have a negligible effect on the stellar distribution in the 
LMC and that general galactic dynamics is the major driver 



of structural change. It is worth noting that our results are 
not necessarily in contradiction with the 'infant mortality' 
scenario, but they show that this mechanism does not ex- 
plain the observed evolution of stellar structure in galaxies. 

Only a small fraction of stars are born in stellar clusters 
which will survive long enough to be able to be identified as 
'clusters' in optical surveys (e.g. Lada & Lada 2003, Gieles 
& Bastian 2008, Bastian 2008). Whether this is due to a high 
fraction of stars born in embedded 'clusters' which quickly 
disrupt or due to a large fraction of stars forming outside 
clusters (or more likely, some combination of these two ex- 
tremes) is still an open question at the moment. 
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APPENDIX A: EFFECTS OF ELLIPTICITY ON 
Q 

In this section we investigate the influence of the elongation 
on the calculated Q parameter. For consistency, we define 
the elongation in the same way as Schmeja & Klessen (2006), 
namely as the ratio of the derived radius when approximat- 
ing the region as a circle to the radius when using the convex 
hull method. Schmeja & Klessen (2006) give the relation be- 
tween the elongation and the major to minor axis ratio of 
an ellipse. 

In order to test the effect we first create 3D spherical 
models of varying fractal dimensions, uniform distribution 
and various centrally concentrated profiles, similar to what 
was done in Cartwright & Whitworth (2004). We then mul- 
tiply one of the three principle axes by a factor and mea- 
sure the Q value and elongation along the three projected 
axes. The results for two fractal models (fractal dimensions 
of 2.0 and 2.5), a uniform distribution and a centrally con- 
centrated distribution with power-law radial index of -2 are 



shown in Fig. IA1I It is clear that increasing the elongation 
decreases the measured Q value. We have also fit each rela- 
tion with a linear slope shown as a dashed line. The error 
bars represent the standard deviation of three realisations 
of the simulations (the fractal models begin with a random 
seed) . 

While the slope of the change in Q as a function of the 
elongation is somewhat dependent on the starting Q value 
(i.e. what type of distribution is used) the overall trend is 
quite similar for all profiles considered. As such, we use the 
average slope of the distributions in the left panel of Fig. I All 
in order to provide a correction function. For this we find 
^intrinsic = Qobs - (~0.208 * elongation) - 0.229. The right 
panel of Fig. lAll shows the resulting values of the simulations 
(shown in the left panel) after this correction has been ap- 
plied. Here the dashed lines represent the input Q value (i.e. 
that determined for the profile where no elongation has been 
applied). While the given correction appears to correct well 
the simulations we note that large errors can be expected 
for systems with high elongations (> 1.5, i.e. major/minor 
axis ratios greater than 2). 



APPENDIX B: FRACTAL MODELS: 
EXTINCTION EFFECTS ON Q 

An additional affect which may influence the measured Q 
parameter is that of differential extinction. This can be real- 
ized conceptually by imagining a uniform stellar population 
and putting a highly sub-structured distribution of gas/dust 
in the line of sight. If the opacity is high enough, this will 
remove whole sections of the stellar distribution from sight, 
creating the appearance of a non-uniform stellar population. 
Here, we we attempt to quantify this affect. 

For the simulations we make 3D fractal distributions of 
point sources with dimension 2.0 following the method pre- 
sented in Goodwin and Whitworth (2004). For each simula- 
tion we create 20,000 gas particles and 1,000 stellar particles 
independently. We then chose a line of sight and counted the 
number of gas particles between the observer and each star 
within some projected tolerance, i.e. we find the surface den- 
sity of gas between the observer and each star. Stars which 
were found to have gas surface densities above a certain 
critical value were removed from the dataset. Finally, we 
calculate Q for the original (non-extincted) and extincted 
datasets and compare their values. Three projections of one 
of the simulations are shown in Fig. IB II where the contours 
represent the integrated column density of gas along the line 
of sight, red circles are stars which have been removed from 
the original dataset (i.e. are found to be extincted), and 
green circles represent stars which were not extincted (i.e. 
the red and green circles together make up the initial stellar 
distribution) . 

The main result is that extinction from clouds tends to 
lower the value of Q, where the change is related to the num- 
ber of stars removed (i.e. the opacity of the clouds). Typical 
changes in the derived Q value, AQ, are lowered by 0.04 to 
0.08 when 20% or 50% of the sources are extincted out of the 
sample, respectively. This makes the distribution appear to 
have a lower fractal dimension than the actual distribution. 
We have tested the effect of varying the number of particles 
used and the input fractal dimensions of the distribution 
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Figure Al. The effect of elongation on Q. Right: To correct for the effect of elongation wc have taken the average slope of the 
simulations (shown in the left panel). Doing this we find Qi n t r i n sic = Qobs ~ (—0.208 * elongation) — 0.229. 



and found that both effects are minimal. Finally, we test 
the effect of culling the stars randomly and find that this 
does not change the value of Q (AQ < 0.015). 

Since an unphysically large number of stars must be 
extincted out of our sample before Q is affected, we conclude 
that extinction does not significantly affect our results. 



APPENDIX C: FRACTAL MODELS: EFFECTS 
OF CONTAMINATION 

Due to a variety of effects, some contamination may be ex- 
pected to enter observational datasets. This could be due 
to field stars in the Galactic halo or young (lower-mass) 
stars which have colours and magnitudes which fall in the 
CMD boxes used to selected older age stars. In order to 
estimate the strength of these effects we have constructed 
the following experiments. We start with a 2D projection 
of a 3D fractal model, chosen to have a similar Q-value as 
our youngest stellar age-box. We then replace a fraction of 
the model stars, whose positions are chosen according to 
the probability function of a background distribution. For 
the background we chose three distributions (all 2D distri- 
butions); a stochastic distribution, and two centrally con- 
centrated distributions, both power-laws, with indices 0.3 
and 1.0. We used the prescription in Cartwright & Whit- 
worth (2004) to construct the background distributions. The 
centrally concentrated distribution with index 0.3 was cho- 
sen as this matches the projected radial profile of the LMC, 
when restricted to the inner three degrees in radius. 

The results of the simulations are shown in Fig. IC 11 
In each case, we see the Q-value increase rapidly with in- 
creasing contamination. In particular, the Q-value in the 



centrally concentrated case (with power-law index of 0.3) 
becomes largely flat for contamination fractions above 0.6. 

This explains why small amounts of fractal substructure 
in the older age boxes (due to young stars which enter our 
selection box) do not significantly affect the results, as a 
large fraction of these stars would need to be present. 

Figure [C2l shows how the two-point correlation function 
analysis is affected by the same contamination. As seen in 
the Q-distribution, the slope and zero-point of the TPCF 
smoothly transform from the value expected for a fractal 
distribution, to that expected for the background or refer- 
ence distributions. 

Thus we conclude that while younger stars (which dis- 
play more substructure) will contribute to all age boxes, this 
effect will not significantly alter our conclusions, as the ex- 
pected contamination is low, their distribution will not sig- 
nificantly change the measured Q distribution or alter the 
TPCF. ' 

This paper has been typeset from a Tj^X/ DTgX file prepared 
by the author. 
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Projected Axis Projected Axis Projected Axis 

Figure Bl. Three projections of a simulation on the effect of extinction on the measured value of Q. Both the gas particles and the 
stars are given a three dimensional fractal distribution of dimension 2.0. The (blue) contours represent the number density distribution 
of gas particles in the simulations (integrated along the line of sight) while the green and red circles represent stars (i.e. together they 
make up the initial sample). The red circles are stars which have been extincted and hence are not 'detected' while the green circles are 
stars which are not extincted enough to be removed from the sample (see text for details). The full sample of stars has a Q value which 
is 0.07 ± 0.03 higher than that of the detected sample. 
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Figure CI. The change in the Q-value of a fractal distribution 
when a given fraction of the stars are replaced by a background 
distribution. Three examples are given for the background distri- 
bution, 1) a stochastic distribution, 2) a centrally concentrated 
power-law with index 0.3 (chosen to match the light profile of the 
LMC in the inner three degrees) and 3) a centrally concentrated 
power-law with index 1.0. 
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Figure C2. The change in the slope and zero-point of the two- 
point correlation function of a fractal distribution when a given 
fraction of the stars are replaced by a background distribution. 
Two examples are given for the background distribution, 1) a 
stochastic distribution, 2) a centrally concentrated power-law 
with index 0.3 (chosen to match the light profile of the LMC 
in the inner three degrees). 
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